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Variability is a key characteristic of young stars. Two major origins may be distinguished: a scaled-up version of the 
magnetic activity seen on main-sequence stars and various processes related to circumstellar disks, accretion and outflows. 
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1 Introduction 

During the first few Myrs of their evolution late-type stars 
pass through several phases that are characterized by the 
amount of circumstellar material and its interaction with 
the forming star. In the earliest protostellar stages the star 
is surrounded by an optically thick envelope (Class O/I ob¬ 
jects). Subsequently the envelope disperses and settles into a 
disk (Class II objects) from which the star may accrete mat¬ 
ter (classical T Tauri star). Finally all circumstellar matter is 
gone and only the star is left, referred to as a Class III ob¬ 
ject or weak-line T Tauri star. The classification with roman 
numbers 0 through III is related to the slope of the spec¬ 
tral energy distribution in the infrared (IR) which is sensi¬ 
tive to the presence of envelopes and disks, while the dif¬ 
ferent T Tauri types represent a spectroscopic classihcation 
measuring accretion through the strength of optical emis¬ 
sion lines; see Feigelson & Montmerle (1999) for details. 
The definition of a ’’young stellar object” (YSO) sometimes 
comprises only the earlier evolutionary phases which are 
dominated by their circumstellar environment, while I in¬ 
clude all the pre-main sequence phases mentioned above in 
this review on YSO variability. 

The variability of YSOs has a variety of origins, last 
but not least due to the complex and diverse environments 
that these objects host. In analogy to the Sun and late-type 
main-sequence stars in general, YSOs display a range of 
magnetic activity phenomena. These comprise coronal ra¬ 
diation from the radio and the X-ray bands, and optical 
and ultraviolet (UV) chromospheric emission lines. Stel¬ 
lar activity goes back to a dynamo mechanism, driven by 
convection and rotation (Noyes et al. 1984). Therefore, the 
faster rotation of young stars is associated with scaled-up 
activity, e.g. 10^ times higher X-ray luminosity is detected 
from YSOs with respect to main-sequence stars (Stelzer & 
Neuhauser 2001; Preibisch & Feigelson 2005). Flares, rep¬ 
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resenting rapid release of magnetic energy through recon¬ 
nection events, are the most prominent activity signatures 
(Sect. Ol. Another ubiquituous type of variability related 
to magnetic activity is rotational modulation produced by 
inhomgeneous surface brightness related to the presence of 
cool star spots (Sect. I2.21 l. while dynamo cycles (Sect 12.31 1 
are observed mostly on old, inactive stars. Next to these 
activity-related variability features, the presence of a cir¬ 
cumstellar disk and mass accretion from the disk can in¬ 
duce variability. The matter infall from the disk to the star 
occurs along magnetic held lines (Koenigl 1991). During 
this process the kinematic energy is transformed into heat, 
leading to the formation of shocks and hot spots upon im¬ 
pact on the stellar surface. The circumstellar environment 
can give rise to variability either due to intrinsic changes in 
mass transfer or disk emission or due to geometric effects as 
accretion streams or structures in the disk rotate in and out 
of the line-of-sight (Sect, [ni l. The most extreme accretion 
events are the sporadic and poorly understood bursts known 
as EXOr and FUOr phenomena which are likely responsible 
for assembling a signihcant amount of the hnal stellar mass 
(Sect l3.21 l. Finally, YSOs are subject to highly dynamic pro¬ 
cesses of mass loss (Sect. [33T l. 

In the remainder of this paper I review recent develop¬ 
ments in variability studies of the above-mentioned origins. 

2 Magnetic activity 
2.1 Flares 

Stellar Hares are the result of magnetic reconnection in the 
upper atmosphere which leads to particle acceleration and 
plasma heating. Footpoints of magnetic loops light up in 
optical/UV emission lines and continuum following elec¬ 
tron bombardment of the lower atmosphere. Intense X-ray 
emission goes along with such events, and is ascribed to the 
Hlling of the coronal loops with heated plasma evaporated 
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off the chromospheric layers along the magnetic field lines 
(‘chromospheric evaporation’, Antonucci et al. 1984). 

The flares observed on YSOs are orders of magnitude 
more energetic than those of main-sequence solar-like stars 
and dMe stars (so-called ’’flare stars”). Typical YSO flare 
energies in the optical broad-bands exceed 10^® erg (Stelzer 
et al. 2003), much larger than the threshold for what is 
defined a ’’super-flare” on a solar-type main-sequence star 
(Maehara et al. 2012). In the X-ray band, similarly high 
flare energies are observed. However, the frequency of large 
flares is small (^ 1 every 10 d for YSO X-ray flares; Stelzer 
et al. 2007a), and sensitivity limits prevent the observation 
of smaller events. Therefore, despite the fact that the ba¬ 
sic picture of the flare physics has been devised decades 
ago (e.g. Cargill & Priest 1983), up to our days simulta¬ 
neous multi-wavelength observations of stellar flares have 
remained extremely rare both for YSOs and for main- 
sequence stars. Exceptions include the first detection of a 
YSO flare at millimeter wavelengths by Bower et al. (2003), 
observed also at centimeter bands and with serendipitous si¬ 
multaneous X-ray coverage. 

If continuous temporal coverage and high enough signal 
for spectral analysis is given hydrodynamical (HD) models 
of flaring plasma loops (Reale et al. 1997) can be adapted 
to the observed time-evolution of X-ray flares, yielding the 
size of the X-ray emitting structures (Getman et al. 2008a). 
Very large loop lengths, multiples of the stellar radius, have 
been found for YSO flares, and it was suspected that these 
events represent magnetic loops that connect from the star 
to the disk in Class II objects (Favata et al. 2005). However, 
after the objects’ evolutionary stage was better constrained 
the evidence for X-ray emitting star-disk loop structures has 
weakened (Getman et al. 2008b; Aarnio et al. 2010). 

The number distributions of flares are a crucial diagnos¬ 
tic for the importance of nano-flares in the coronal heat¬ 
ing process. They have been extensively studied for the 
Sun across a wide wavelength range from the UV to the 
hard X-ray band (see e.g. references in Aschwanden 2014). 
For YSOs the X-ray flare energies follow a power-law with 
spectral index Ri 2 (Stelzer et al. 2007a; Albacete Colombo 
et al. 2007; Caramazza et al. 2007), suggesting sufficient 
energy in small events to explain the existence of the X-ray 
corona through nano-flaring. However, the YSO X-ray flare 
energy distributions show a turn-over for energies below 
~ 10^® erg/s due to incompleteness (Fig[T]), such that the 
nano-flare heating hypothesis relies on extrapolation over 
many orders of magnitude. Analogous to the X-ray band, 
flare energy number distributions have been studied for a 
small number of late-type main-sequence stars in the ex¬ 
treme UV (Audard et al. 2000) and optical bands (Hawley 
et al. 2014) but equivalent studies at those wavelengths are 
not yet available for YSOs. 

The importance of studying X-ray flares goes beyond 
the scope of understanding coronal physics. The UV and 
X-ray emission of YSOs is a major heating and ionization 
agent of their circumstellar disks (e.g., Gorti & Hollenbach 



log Ep [erg] 


Fig.l Cumulative number distributions of X-ray flare en¬ 
ergies for YSOs in the Orion Nebula Cluster (ONC) ob¬ 
served during the Chandra Orion Ultradeep Project (COUP) 
and in the Taurus Molecular Clouds (TMC) observed during 
the XMM-Newton Extended Survey of the Taurus Molecular 
Clouds (XEST); from Stelzer et al. (2007a). 

2009; Ercolano et al. 2008). Higher-energy radiation pene¬ 
trates deeper and, therefore, the effects on the disk crucially 
depend on the intensity and spectral signature of the stellar 
radiation (Ercolano & Glassgold 2013). 

2.2 Rotational modulation 

Periodic optical variability of stars on time-scales of days is 
explained as modulation of the surface-averaged brightness 
when cool and dark star spots move across the line-of-sight. 
Next to flares, these are the most frequent and well-studied 
stellar variability phenomenon and a means to measure the 
stellar rotation rates (see e.g. Herbst et al. 2000). However, 
in Class II YSOs the optical emission and its variability is 
often dominated by processes related to their circumstellar 
material from the star formation process (see Sect.O. 

In the X-ray band, rotational modulation associated with 
structured coronae can be expected but is usually masked 
by much stronger irregular variability related to flaring. A 
clear periodic X-ray signal persisting over several years has 
been observed on the Class I protostar VI647 Ori, an EX Or 
type variable (see Sect. l3.21 i. and the phase-folded lightcurve 
was modelled with emission from accretion hot spots (Ham- 
aguchi et al. 2012). In fact, hot, bright spots representing 
the impact points of accretion streams on the stellar surface 
may in principle produce similar rotational signals to dark 
magnetic spots (see Sect l3.1l) . 

An nearly uninterrupted 13 d-long Chandra observation 
- COUP, introduced in Fig. [T] - has enabled a search for 
X-ray variability on rotational timescales of YSOs through 
periodogram analysis. X-ray periods were found in 10 % 
of the studied sample (Flaccomio et al. 2005). Statisti¬ 
cal methods examining changes of the amplitude of X-ray 
lightcurves throughout the phases of the known periodic op¬ 
tical variability pattern have also been used to establish the 
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presence of rotational modulation in the X-ray band (Flac- 
comio et al. 2012). While amplitude peaks at optical phases 
(p ~ 0.5 and 1.5 were observed for all types of YSOs, they 
were most pronounced for Class II objects; see Sect. l3.1l for 
details. 

2.3 Dynamo cycles 

Dynamo cycles, similar to the Sun’s 11-yr cycle, have been 
observed on ~ 60% of solar-type main-sequence stars 
through variability of the Ca II H&K index (Baliunas et al. 
1995). For a much smaller number of late-type stars, dy¬ 
namo cycles have been identified through photometric mon¬ 
itoring (e.g. Olah et al. 2009), including some young main- 
sequence stars. However, only one dynamo cycle is known 
on a YSO so far, the weak-line TTauri star V410Tau in 
the ^ 2Myr-old Taurus molecular clouds. The cycle of 
V410Tau was detected in the V band and has a relatively 
short period of 5.4 yrs (Stelzer et al. 2003, Olah et al. 2009). 

While broad-band photometry traces star spots, calcium 
emission is sensitive to chromospheric structures such as 
plages. The existing data may suggest that dynamo cycles 
on younger stars are predominantly visible through spot- 
variations, while cycles on more evolved stars are better 
probed through chromospheric tracers. However, strong bi¬ 
ases - in particular the limited amount of time invested in 
the required long-term monitoring observations - may af¬ 
fect our picture. For the same reason, unambiguous long¬ 
term periodic variability in the X-ray band has been iden¬ 
tified on only four stars so far. The youngest of these is 
ss 500 Myr old (Sanz-Forcada et al. 2013), clearly beyond 
the YSO stage. A possible dependence of the cycle period 
on age, i.e. on the overall strength of activity, needs further 
investigation. 

3 Circumstellar environment 

3.1 Accretion and disks 

The prime diagnostics for YSO disks and accretion are the 
IR and optical wavelengths. While disks are bright in the 
IR, the accretion process gives rise to optical and IR emis¬ 
sion lines. Searching for the origin of the radiation emit¬ 
ted from star-disk systems is best pursued in the time do¬ 
main, as different variability patterns and timescales are 
expected for rotational modulation and for intrinsic vari¬ 
ations. Magneto-HD simulations predict for the stable ac¬ 
cretion mode periodic broad-band lightcurves and emission 
line characteristics resulting from the changing visibility of 
accreting matter in the magnetic funnel streams throughout 
a stellar rotation (Romanova et al. 2008; Kurosawa et al. 
2008). These variations occur on timescales of the YSO ro¬ 
tation period. Quasi-periodic variability on time-scales of 
several rotation cycles is expected in the case of structural 
changes of the magnetosphere as result of star-disk recon¬ 
nection (Goodson et al. 1997). Intrinsic variations may be 



Fig. 2 Near-IR color-color diagram for the TWA 30 AB 
binary with data from Looper et al. (2010) and new observa¬ 
tions marked with filled symbols; Principe et al., ApJ subm. 

caused by a variety of processes including unstable accre¬ 
tion through the Rayleigh-Taylor instability. The timescales 
for these processes are much shorter than a rotation cycle 
and they produce stochastic variability patterns in photo¬ 
metric lightcurves and emission line profiles (Romanova et 
al. 2008; Kurosawa et al. 2013). Both photometric and spec¬ 
troscopic monitoring are, therefore, powerful tools to distin¬ 
guish different variability mechanisms in YSOs. 

Wolk et al. (2013) showed how a combined analysis of 
near-IR magnitude and color changes and their time-scales 
helps to unveil the origin of the variability: Variations pro¬ 
duced by both hot and cool spots are periodic on timescales 
of the stellar rotation cycle but show opposite color trends; 
extinction variations are non-periodic or associated with the 
disk rotation and they tend to induce dips in the lightcurve; 
variations of the mass accretion rate induce strong bright¬ 
ness changes. A case of very different near-IR variability in 
two stars of a young binary is TWA 30 AB (Looper et al. 
2010). While the variations observed in TWA30A follow 
the extinction vector suggesting obscuration effects by the 
disk, the near-IR color changes of TWA 30B occur along the 
locus populated by classical T Tauri stars indicating varying 
amount of dust emission from the disk. However, new ob¬ 
servations challenge this clear distinction (see Fig.lHl. 

The advent of space-borne optical and IR monitoring 
instruments has boosted research of photometric YSO vari¬ 
ability associated with disks and accretion. Weeks or even 
months of uninterrupted data have become available from 
space missions such as Spitzer and CoRoT (e.g., Morales- 
Calderon et al. 2011; Cody et al. 2014). The morphology 
of irregular variations allows e.g. to distinguish stellar flares 
from accretion bursts (Stauffer et al. 2014). Performing such 
observations for a significant number of YSOs - as can 
readily be achieved with the large field-of-views of imag¬ 
ing instruments given the clustered nature of YSOs - yields 
information on the relative frequencies of different types of 
variability patterns (Cody et al. 2014). 

Variable mass accretion, whether intrinsic or caused by 
geometric effects, has been cited as a possible cause of 
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the large observed spread in the relation between mass ac¬ 
cretion rate (Macc) and stellar mass (e.g. Mohanty et al. 
2005). However, recent results on a large sample of accret¬ 
ing YSOs indicate that the spread of Macc for given mass 
is signihcantly larger than what can be accounted for by 
variability (Venuti et al. 2014). The accretion rates deter¬ 
mined in this study have been obtained from photometric 
measurements of UV excess emission shown by Gullbring 
et al. (1998) to be correlated with Macc- Alternative diag¬ 
nostics for Macc are optical and near-lR emission lines. The 
Ha line is usually the brightest emission feature in YSO 
spectra, and therefore the most widely used probe of ac¬ 
cretion despite the fact that its profile has long been rec¬ 
ognized to comprise contributions from outflows (Calvet 
1997). This is the likely reason why accretion rate esti¬ 
mates from Ha emission yield vastly different Macc values 
than other prominent accretion signatures, such as higher-n 
Balmer lines or calcium (Ca II H&K and Ca II IRT) and he¬ 
lium (e.g.. He IA587, He lA1083) lines (Alcala et al. 2014). 

For a given star, the individual emission lines can show 
different variability patterns (e.g. Donati et al. 2008; Stelzer 
et al. 2014) suggesting different formation regions. In a 
long-term, low-cadence spectroscopic monitoring program 
of a small sample of YSOs the amplitude of emission line 
variability reached a maximum after 25 d, indicating that the 
dominant timescales for accretion variability is < 2 — 3 
weeks (Figl2 Costigan et al. 2012). On a shorter observ¬ 
ing cadence, both slow line profile variations suggestive of 
rotational modulation and rapid events possibly related to 
disk or magnetosphere instabilities were observed (Costi¬ 
gan et al. 2014). Similarly, dedicated Ha monitoring of 
the young benchmark brown dwarf 2MASSW J1207334- 
393254 showed both smooth variations on a timescale cor¬ 
responding to the rotation period and fast variations imply¬ 
ing drastic changes of the mass accretion rate (Stelzer et al. 
2007b; Scholz et al. 2005). The similarity of the variabil¬ 
ity pattern with stellar-mass accreting YSOs underlines that 
substellar objects undergo a T Tauri phase. 

Multi-wavelength monitoring combining photometric 
and spectroscopic observations has the potential to connect 
variations related to accretion (emission lines) to those re¬ 
lated to disk emission (IR brightness). Such campaigns are 
technically challenging, have therefore been limited to few 
objects and have so far not provided clear results (Faesi et 
al. 2012) except for individual cases. One of the best test- 
cases for YSO variability is AA Tau, a classical T Tauri star, 
which thanks to its known, high inclination (i ^ 75°) is 
ideally suited to demonstrate the complexity of the YSO cir- 
cumstellar environment. Among the characteristic features 
of AA Tau are (i) rotationally modulated veiling variations 
indicating the presence of hot, accretion spots, (ii) time- 
delays between emission line peaks and veiling variations, 
(iii) photometric brightness variations interpreted as occul- 
tations by a disk warp, and (iv) a recent fading of the op¬ 
tical brightness without associated changes of its Macc and 
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Fig. 3 Average spread in accretion rate for 8 Class II 
YSOs in the Chal star forming region over various 
timescales showing little increase beyond the shortest 
timescale probed (8...25d); from Costigan et al. (2012). 

attributed to enhanced dust extinction from a density pertur¬ 
bation in the disk (Bouvier et al. 2003; Bouvier et al. 2013). 

There are manifold interactions between the YSO X- 
ray emission resulting from magnetic activity on the stel¬ 
lar surface and the circumstellar environment. First, as 
mentioned in Sect. 12.21 differences have been observed in 
the X-ray variability pattern of YSOs with and without 
disks/accretion. Flaccomio et al. (2012) ascribed them to 
time-variable absorption of the coronal X-ray emission by 
structures in the circumstellar material, such as disk warps 
or accretion streams. The same phenomenon may also ex¬ 
plain the on average factor-two deficiency of the X-ray lu¬ 
minosity of Class II objects with respect to YSOs without 
disks (Gregory et al. 2007). 

While the presence of a YSO disk may affect the visibil¬ 
ity of the stellar corona as seen above, an influence of X-ray 
emission on disk structure and accretion has been advocated 
in several theoretical works. Ke et al. (2012) predicted the 
enhanced ionization of disk grains resulting from the strong 
X-ray flux during YSO coronal flares to produce changes 
in the height of the inner disk rim and associated IR vari¬ 
ability. In a different scenario, flares in the star-disk mag¬ 
netosphere provoke a density perturbation in the disk which 
gives rise to a mass accretion event (Orlando et al. 2011). 
No convincing correlation between X-ray and IR variabil¬ 
ity in YSOs has been observed so far (Flaherty et al. 2014), 
and the relation between X-ray properties and Macc is un¬ 
explored to date. However, the typical low cadence of X-ray 
observations may have prevented the detection of joint vari¬ 
ations with optical/IR light on short timescales. An excep¬ 
tion is the case of TWHya on which Dupree et al. (2012) 
have identified episodes of enhanced mass accretion in soft 
X-ray emission produced in the accretion shoclfl The ob¬ 
served soft X-ray variability of TW Hya was contemporane- 

* YSO X-ray emission at energies < 0.5 keV may be produced in ac¬ 
cretion shocks which can be heated by infalling matter up to ~ 3 — 4 MK. 
The tendency of Class II objects to show an excess of soft X-rays with re¬ 
spect to Class III objects was understood as evidence for accretion - rather 
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ous with profile variations in Balmer lines indicating down¬ 
flow towards the star, and they were followed by an increase 
in veiling with a few-hours delay compatible with the ex¬ 
pected inflow timescale. Moreover, the increase in veiling 
was correlated with the harder (> IkeV) coronal X-ray 
flux, supporting the picture of an ’’accretion-fed corona” 
where the stellar X-ray emission is enhanced as a result of 
heating by the accreting matter (Brickhouse et al. 2010). 

3.2 FUOr and EXOr events 

A small number of YSOs has been observed to exhibit large 
brightness variations (up to 4 — 5 mag). These objects are 
classified according to the timescales of the outbursts as 
FU Ors (duration of decades; Hartmann & Kenyon 1985) or 
as EX Ors (duration of months and recurrent; Herbig 1989). 
Both types of outbursts are ascribed to drastic accretion rate 
changes. A major part of the stellar mass may be assembled 
during such outbursts but the physical mechanisms trigger¬ 
ing the outburst are still debated; see Audard et al. (2014) 
for a recent review. 

Contradictory results have been obtained from the study 
of YSO X-ray emission and its variability during accretion 
outbursts. While the X-ray lightcurves of EX Ors tend to 
be correlated with the optical/lR flux, no clear picture has 
emerged yet for the characteristics and origin of the X-ray 
emission, and only a handful of objects have been moni¬ 
tored in X-rays. In most of these, the temperature of the X- 
ray emitting gas was found to be too high for being gener¬ 
ated in accretion shocks, but see Teets et al. (2012) for a dif¬ 
ferent example. An alternative interpretation for the soften¬ 
ing of the X-ray emission observed during some optical/IR 
outbursts is a change in the structure of the X-ray emitting 
magnetosphere caused by the accretion event (Audard et al. 
2005). 

3.3 Jets and outflows 

Mass ejection is an essential ingredient of pre-main se¬ 
quence evolution and might be responsible for carrying 
away the excess angular momentum of the accreted matter. 
The highly collimated, fast moving (~ 100 — 1000 km/s) 
inner parts are referred to as ’’jets” and the wide-angle 
slowly moving ('^ 10 km/s) molecular ejecta as ’’outflows”; 
see Erank et al. (2015) for a review. YSO jets evolve on 
timescales of few years and are, thus, easily accessible to 
variability studies. Since the first proper motion measure¬ 
ments of Herbig & Jones (1981), the bow-shocks associ¬ 
ated with YSO jets - so-called Herbig-Haro objects - have 
been imaged for YSOs in a wide range of evolutionary 
stages (Class 0 to Class II). In particular, narrow-band imag¬ 
ing with the Hubble Space Telescope thanks to its high spa¬ 
tial resolution has enabled detailed studies of the morphol¬ 
ogy and dynamics of YSO outflows (see e.g. Hartigan et al. 

than the corona - playing a major role in this X-ray component (Giidel & 
Telleschi 2007). 



Fig. 4 Chandra/ACIS images for the YSO ZCMa in the 
0.2 — IkeV band obtained in Dec 2003 (left) and in Dec 
2008 (right). The detected X-ray sources are marked with 
black circles. The unresolved binary components of ZCMa 
are indicated with white x-shaped symbols. The position of 
the X-ray source detected in the 2008 image to the SW of 
Z CMa, labeled ‘J’, is overlaid in the 2003 image as red cir¬ 
cle. In 2008, the soft X-ray image shows elongated emission 
along the optical jet axis. (Eigure adapted from Stelzer et al. 
2009). 

2011). Next to the propagation of individual knots, variabil¬ 
ity due to jet precession or orbital motion in binaries can 
be tracked through optical and IR imaging by studying the 
displacements of different knots from a linear axis pointing 
towards the star (e.g. Cunningham et al. 2009). These vari¬ 
ations occur on secular timescales and can not be traced by 
multi-epoch observations. 

Roughly a dozen YSO jets have been detected in X-rays 
so far (see summary in Bonito et al. 2007). Variations of the 
source positions and strength in X-ray images taken a few 
years apart were observed in a few of these (Eavata et al. 
2006; Stelzer et al. 2009 and Fig.|4l Giidel et al. 2011). For 
typical shock velocities on the order of 500 km/s, jet proper 
motion can be detected in the most nearby YSOs (~ 100 pc) 
in less than a decade. Yet it is unclear to date if the X-ray 
knots represent propagated plasma ejected and heated sev¬ 
eral years prior to the observation or transient emission from 
locally heated physical structures (Schneider et al. 2011). 
Short predicted cooling times (Guenther et al. 2009) favor 
the latter interpretation. On the other hand, the outward de¬ 
crease of density results in a decrease of heating efficiency 
such that X-ray jets are expected to be most luminous close 
to the star. In fact, the X-ray jets imaged so far were found at 
separations < 5" from the central driving source, resolvable 
at present only with Chandra. 

The observational evidence of both stationary and mov¬ 
ing components in X-ray jets has been addressed in two- 
dimensional HD simulations in which the X-ray emission 
from a supersonic jet interacting with its ambient medium 
is synthesized. Bonito et al. (2010) showed that a pulsed jet 
with different values for crucial jet parameters, such as den¬ 
sity and ejection rate, can describe the X-ray luminosities 
and temperatures of some of the observed YSO jets but it 
fails to explain the low X-ray temperatures of the knots with 
the largest separation from the star (> 1000 AU). The sta¬ 
tionary (non-variable) X-ray component seen in some jets 
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can be reproduced if a nozzle at the base of the jet is in¬ 
troduced into the model (Bonito et al. 2011). Currently our 
understanding of the high-energy component of YSO jets 
is limited mainly by the large uncertainties of the observed 
X-ray parameters imposed by sensitivity constraints. 


References 

Aarnio, A.N., Stassun, K. G. & Matt, S.P., 2010: ApJ 717, 93 
Albacete Colombo, R, Caramazza, M., Flaccomio, E., Micela, G. 

& Sciortino, S., 2007: A&A 474, 495 
Alcala, J.M., Natta, A., Manara, C.F., et al.: 2014, A&A 561A, 2 
Antonucci, E., Gabriel, A. H., Dennis, B. R.: 1984, ApJ 287, 917 
Aschwanden, M.J., 2014: ApJ 782, 54 

Audard, M., Abraham, R, Dunham, M.M., et al.: 2014, Protostars 
& Planets VI, University of Arizona Press, p.387 
Audard, M., Gudel, M., .Skinner, S.L., et al.: 2005, ApJ 635, L81 
Audard M., Gudel M., Drake J.J. & Kashyap V.L., 2000: ApJ 541, 
396 

Baliunas, S.L., Donahue, R.A., Soon, W.H., Home, J.H., Frazer, 
J., Woordard-Eklund L., et al.: 1995, ApJ 438, 269 
Bonito, R., Orlando, S., Miceli, M.: 2011, ApJ 737, 54 
Bonito, R., Orlando, S., Miceli, M., et al.: 2010, A&A 517A, 68 
Bonito, R., Orlando, S., Peres, G., Favata, F. & Rosner, R.: 2007, 
A&A 462, 645 

Bouvier, J., Grankin, K., Ellerbroek, L. E., Bouy, H. & Barrado, 
D.:2013, A&A557A, 77 

Bouvier, J., Grankin, K.N., Alencar, S.H.P, et al.: 2003, A&A 409, 
169 

Bower, G.C., Plambeck, R.L., Bolatto, A., et al.: 2003, ApJ 598, 
1140 

Brickhouse, N.S., Cranmer, S.R., Dupree, A.K., Luna, G.J.M. & 
Wolk, S.: 2010, ApJ 710, 1853 
Calvet, N.: 1997, lAU Symposium 182, 417 
Caramazza, M., Flaccomio, E., Micela, G., Reale, R, Wolk, S. J. 

& Eeigelson, E. D., 2007: A&A 471, 645 
Cargill, PJ. & Priest, E.R.: 1983, ApJ 266, 383 
Cody, A.M., Stauffer, J., Baglin, A., et al.: 2014, AJ 147, 82 
Costigan, G., Vink, J. S., Scholz, A., Ray, T. & Testi, L.: 2014, 
MNRAS 440, 3444 

Costigan, G., Scholz, A., Stelzer, B., Ray, T., Vink, J.S. & Mo- 
hanty, S.: 2012, MNRAS 427, 1344 
Cunningham, N.J., Moeckel, N. & Bally, J.: 2009, ApJ 692, 943 
Donati, J.-R, Jardine, M.M., Gregory, S.G., et al.: 2008, MNRAS 
386, 1234 

Dupree, A.K., Brickhouse, N.S., Cranmer, S.R., et al.: 2012, ApJ 
750, 73 

Ercolano, B. & Glassgold, A.E.: 2013, MNRAS 436, 3446 
Ercolano, B., Drake, J.J., Raymond, J.C. & Clarke, C.C.: 2008, 
ApJ 688, 398 

Faesi, C.M., Covey, K.R., Gutermuth, R., et al.: 2012, PASP 124, 
1137 

Favata, R, Bonito, R., Micela, G., et al.: 2006, A&A 450, L17 
Favata, R, Flaccomio, E., Reale, R, et al.: 2005, ApJS 160, 469 
Eeigelson, E.D. & Montmerle, T.: 1999, ARA&A 37, 363 
Flaccomio, E., Micela, G. & Sciortino, S.: 2012, A&A 548A, 85 
Flaccomio, E., Micela, G. & Sciortino, S., et al.: 2005, ApJS 160, 
450 

Frank, A., Ray, T., Cabrit, S., et al.: 2014, Protostars & Planets VI, 
University of Arizona Press, p.451 
Flaherty, K.M., Muzerolle, J., Wolk, S.J., et al.: 2014, ApJ 793, 2 


Getman, K.V, Eeigelson E.D., Broos P.S., Micela G. & Garmire 
G.P.: 2008a, ApJ 688,418 

Getman, K.V, Eeigelson ,E.D., Micela, G., Jardine, M.M., Gre¬ 
gory, S.G., & Garmire, G.P: 2008b, ApJ 688, 437 
Goodson, A. P, Winglee, R. M. & Bdhm, K.-H.: 1997, ApJ 489, 
199 

Gorti, U. & Hollenbach, D.: 2009, ApJ 690, 1539 
Gregory, S.G., Wood, K. & Jardine, M.: 2007, MNRAS 379, L35 
Gudel, M., Audard, M., Bacciotti, R, et al.: 2011, ASPC 448, 617 
Gudel, M. & Telleschi, A.: 2007, A&A 474, L25 
Gullbring, E., Hartmann, L., Briceno, C. & Calvet, N.: 1998, ApJ 
492, 323 

Gunther, H.M., Matt, S.P & Li, Z.-Y.: 2009, A&A 493, 579 
Hamaguchi, K., Grosso, N., Kastner, J.H., et al.: 2012, ApJ 754, 
32 

Hartigan, P, Frank, A., Foster, J.M., et al.: 2011, ApJ 736, 29 
Hartmann, L. & Kenyon, S.J.: 1985, ApJ 299, 462 
Hawley, S.L., Davenport, J.R.A., Kowalski, A.R, Wisniewski, J.P, 
et al.: 2014, ApJ 797, 121 

Herbig, G.H.: 1989, ESO Workshop on Low Mass Star Formation 
and Pre-Main Sequence Objects, p.233 
Herbig, G.H. & Jones, B.R: 1981, AJ 86, 1232 
Herbst, W., Maley, J.A. & Williams, E.C.: 2000, AJ 120, 349 
Ke, T.T., Huang, H. & Lin, D.N.C.: 2012, ApJ 745, 60 
Koenigl, A.: 1991, ApJL 370, L39 

Kurosawa, R. & Romanova, M.M.: 2013, MNRAS 431, 2673 
Kurosawa, R., Romanova, M.M. & Harries, T.J.: 2008, MNRAS 
385, 1931 

Looper, D.L., Bochanski, J.J., Burgasser, A.J., et al.: 2010, AJ 140, 
1486 

Maehara, H., Shibayama, T., Notsu, S., et al.: 2012, Nature 485, 
478 

Mohanty, S., Jayawardhana, R. & Basri, G.: 2005, ApJ 626, 498 
Morales-Calderon M., Stauffer, J.R., Hillenbrand, L.A., et al.: 
2011, ApJ 733, 50 

Noyes, R.W., Hartmann, L.W., Baliunas, S.L., Duncan, D.K. & 
Vaughan, A.H.: 1984, ApJ 279, 763 
Olah, K., Kollath, Z., Granzer T., Strassmeier K.G., Lanza, A.R, 
Jarvinen, S., et al.: 2009, A&A 501, 703 
Orlando, S., Reale, R, Peres, Go. & Mignone, A.: 2011, MNRAS 
415, 3380 

Preibisch, T. & Eeigelson, E.D.: 2005, ApJS 160, 390 
Principe, D.A., Kastner, J.H., Sacco, G., et al.: 2015, in prep. 
Reale, R, Betta, R., Peres, G., Serio, S. & McTiernan, J.: 1997, 
A&A 325, 782 

Romanova, M.M., Kulkarni, A.K. & Lovelace R.V.E.: 2008, ApJ, 
673L, 171 

Sanz-Forcada, J., Stelzer, B. & Metcalfe, T.S.: 2013, A&A 553, L6 
Scholz, A., Jayawardhana, R. & Brandeker, A.: 2005, ApJL 629, 
41 

Schneider, PC., Gunther, H.M. & Schmitt, J.H.M.M.: 2011, A&A 
530A, 123 

Stauffer, J., Cody, A. M., Baglin, A., et al.: 2014, AJ 147, 83 
Stelzer, B., Alcala, J.M., Whelan, E. & Scholz, A.: 2014, EPJ Web 
of Conferences, Vol.64, id.08005 
Stelzer, B., Hubrig, S., Orlando, S., et al.: 2009, A&A 499, 529 
Stelzer, B., Flaccomio, E., Briggs, K., Micela, G., Scelsi, L., et al.: 
2007a, A&A 468, 463 

Stelzer, B., Scholz, A. & Jayawardhana, R.: 2007b, ApJ 671, 842 
Stelzer, B., Fernandez, M., Costa, V. M., Gameiro, J. R, Grankin, 
K., etal.: 2003, A&A 411, 517 
Stelzer, B. & Neuhauser, R.: 2001, A&A 377, 538 


© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Wdnhdm 


www.an-journal.org 




Astron. Nachr. / AN (2006) 


795 


Teets, W.K., Weintraub, D.A., Raster, J.H., et al.: 2012, ApJ 760, 
89 

Venuti, L., Bouvier, J., Flaccomio, E., et al.: 2014, A&A 570, A82 
Wolk, S.J., Rice, T.S., Aspin, C.: 2013, ApJ 773, 145 


www.an-journal.org 


© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 




